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1.  Introduction 

Hepatitis  B  virus  infection  is  fifty  to  one  hundred  times  more  transmissible  than  HIV  [1],  The 
management  of  HIV-HBV  co-infection  is  difficult  by  the  use  of  drugs  with  action  against  both 
viruses,  the  risk  of  flares  and  hepatic  decompensation  with  immune  reconstitution,  and  the 
increasing  incidence  of  antiviral  resistance.  Due  to  shared  modes  of  transmission,  co-infection  is 
general  and  an  estimated  4  million  people  worldwide  are  co-infected  with  HBV-HIV.  Some  drugs 
used  to  treat  HIV  are  poisonous  to  the  liver,  which  may  amount  to  more  damage  from  the  hepatitis 
B  infection  [7].  Individuals  co-infected  with  both  hepatitis  B  and  HIV  are  14  to  17  times  more 
likely  to  die  than  those  with  hepatitis  B  alone  [6], 

Hepatitis  B  is  a  dynamic  disease  and  it  is  imperative  to  understand  its  virology  and  natural  history 
in  other  to  reduce  complications  and  limit  the  progression  of  disease.  The  use  of  drugs  with 
activity  against  both  viruses  (HIV-HBV  co-infection),  the  risk  of  flares  and  hepatic 
decompensation  with  immune  reconstitution,  and  the  increasing  prevalence  of  antiviral  resistance 
complicate  the  management  of  HIV  -  HBV  co  -  infection.  Forty  million  people  around  the  world 
are  infected  with  HIV.  To  make  things  worse,  certain  medicines  used  to  treat  HIV  are  toxic  to  the 
liver  that  may  already  be  damaged  by  the  infection  with  hepatitis  B.  In  this  study,  a  realistic 
mathematical  model  for  HBV  -  HIV  co  -  infection  is  formulated  and  analyzed,  incorporating  the 
key  epidemiological  and  biological  characteristics  of  each  of  the  two  diseases.  This  study's  main 
contribution  is  to  conduct  an  optimal  test  of  the  resulting  model. 
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2.  Mathematical  Model  Formulation 

The  model  is  constructed  based  on  the  HBV/HIV  transmission  framework  and  the  model  proposed 
in  [6],  The  population  is  divided  into  thirteen  classes  based  on  epidemiological  status;  a  co- 
infections  models  of  susceptible  individuals  of  both  diseases  (S  (t)).  Infected  individuals  in  the 
asymptomatic  stage  of  HIV  infection  (Hi  (t)).  HIV-infected  individuals  that  exhibit  clinical 
symptoms  of  AIDS  (H2  (t)).  Dually-infected  individuals  with  HBV  acute  infection,  in  the 
asymptomatic  acute  stage  of  HIV  infection  (H11  (t)).  Dually-infected  individuals  with  HBV  acute 
infection,  displaying  symptoms  (symptomatic)  of  AIDS  (H21  (t)).  Dually-infected  individuals  with 
HBV  chronic  infection,  in  the  asymptomatic  stage  of  HIV  infection  (Hie  (t)).  Dually-infected 
individuals  with  HBV  chronic  infection,  displaying  symptoms  of  AIDS  (H2C  (t)).  HBV  recovered 
individuals  with  protective  immunity  in  the  asymptomatic  stage  of  HIV  infection  (Hir  (t)).  HBV 
recovered  individuals  with  protective  immunity  displaying  symptoms  of  AIDS  (H2R  (t)).  HBV 
vaccinated  individuals  (V  (t)).  Individuals  with  HBV  acute  infection  (Ib  (t)).  HBV  chronic  carriers 
(Cb  (t))  and  HBV  recovered  with  protective  immunity  (Rb  (t)). 


^  l  l 


Figure  1:  Schematic  representation  of  interactions  of  Co-infection,  HBV-HIV  transmission 
^  =  7t  +  ii(o{\-vCb-v1Hic-v2H2C)  +  tV-{\-8){1h+1b)S-{^  +  8)S 

=  -6))  +  SS  +  crIB-(AH  +  ttbAb)V  -(a,+t)V 

JTJ 

-^  =  (v  +  (l-S)s)AH-<p1(l-a)ABHl-(Ju0+a)H1, 
dH 

—  aH{  —  (p2ABH2  —  +  dH  )  H2  , 


(1) 

(2) 

(3) 
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dHu 

dt 

dH  2  j 
dt 

dHlc 

dt 

dH2  c 
dt 

dHi 

dt 

dH2R 

dt 


=  <p2XBH2  +  ccxHXI  —  (pi0  +  dXH  +  (j>2 )  H21  , 

—  pia)vxIix c  +  qx(f>x  (l —  ccx ) Hl  j  +y/ 2  (\  —  y^}  A,HCB  —  +  d4H  +  <22  +  ^  (l —  oc2 ))  Hlc 

—  /jcov2H2  c  +  q2tj)2H2  I  +  ot2Hlc  —  (/j0  4~  d2H  +  02 )  H2C  , 

—  ^1  (l  —  ^1  )  (^  _  *?1  )  ^1./  ^1  (l  —  ^2  )  ^1,C  ~  (Mo  ^3  )  ^IR  ’ 

—  (l  —  ^2  )^2/  @2^20  ^3^1R  ~  (M)  ^3H  )^2R  ’ 


(5) 

(6) 

(7) 

(8) 

(9) 

(10) 


d[s_ 

dt 

dCH 

dt 


(kbV  +  {\-5)s)XB-y/(l-yl)XHlB-{^  +  c>+yl)lB 


ficovC H  +  cjytI  B  \j/ 2  (l  y2)/^H^-'B  ( /Ai  d B  ^y2) CB 


dR 

——  —  (l  —  q}yxI B  +  y2CB  —  i//2AhRb  —  £iQRB 
dt 


(ID 

(12) 

(13) 


Where, 

-  _  A  (/b  +  Vi  {H\i  +  Vih^ii )  7  (CB  +  t]lcHlc  +  t]2cH2C  )) 


Pu  {H i  +  £H2  +  ( //,/  +  slhH2I )  +  £\  (fflc  +  t:lhIllc  )  +  £3  (771S  +  £yj^2i; )) 


'bH  ~ 

Table  1:  T1 

N 

he  parameters  used 

Parameter 

Interpretation 

Value 

Reference 

M 

Birth  rate  in  the  population 

0.0121 

[5] 

Mo 

Natural  mortality  rate 

0.00693 

[5] 

Mi 

HBV  related  mortality  rate 

0.007 

[4] 

CO 

Proportion  of  births  without  successful  vaccination 

0-100  percent 

[2] 

0 

Vaccination  rate  of  acute  HBV  infectious  class  Ib 

Hypothetical 

PB 

Effective  contact  rate  for  HBV 

0.4 

[5] 

Pu 

Effective  contact  rate  for  HIV 

0.03 

[5] 

K 

The  recruitment  rate  of  susceptible  individuals  in  the  population 

500 

[5] 

kb 

Rate  of  waning  of  vaccine 

0.1 

[5] 

T 

The  rate  of  relapse  of  vaccine  induced  immunity 

Hypothetical 

(px  and  (p2 

Probabilities  of  acquiring  HBV  infection  of  individuals  in  HiandH2 

0.3,  0.5 

[5] 
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VvVi 
and  i//3 

Probabilities  of  acquiring  HIV  infection 
of  individuals  in  Ib  ,  Cb  and  Rb 

0.4,  0.5,  0.1 

[5] 

and  v2 

The  proportion  of  perinatal  infected  bom  by  carrier  mothers  in  the  Cb, 
Hie  and  H2C  classes 

0.11,0.12,  0.13 

[5] 

n,  fa 

and  (j)2 

Rate  of  moving  from  Ib,  Hii  and  H21  classes  to  Cb,  Hie  and  H2C 
classes 

4,  2,3 

[5] 

q,  q2 
and  q3 

The  average  probability  of  moving  from  Ib,  Hu  and  H21  classes  fail  to 
clear  an  acute  state  and  develops  a  carrier  state  Cb,  Hie  and  H2C 
classes 

0.885,  0.9,  0.95 

[5] 

y 2  ’  $ 
and  d2 

Rate  of  moving  from  Cb,  Hie  and  H2C 

Classes  for  Rb,  Hir  and  H2R  classes 

0.025,  0.05,  0.06 

[5] 

dB 

HBV  induced  death  rate  in  Cb  class 

0.002 

[5] 

dH  ?  dlH , 

^2  H  ’  ^3  H 

and  dAH 

HIV  induced  death  rate  of  individuals  in 

H2,  H11,  H21,  Hie,  H2C  and  H2R  classes 

0.01,0.03,0.05, 

0.02,  0.06 

[5] 

8 

The  vaccination  rate 

0.75 

[5] 

Th’Thh’tb 

qlcandq2C 

modification  parameters  of  the  relative  infectiousness  of  individual  in 
the  Cb,  Hu,  H21,  Hie,  H2c 

1.1,  1.2,  0.16,  1.3, 

1.4 

[5] 

£V  £\ h’ 

£1  ’  £2h  ’^3  ’ 

£3h 

modification  parameters  of  the  relative  infectiousness  of  individual  in 
the  H11,  H21,  H2,  Hie,  H2C,  Hir  and,  H2R 

1.2,  1.25,  1.1,  1.5, 
1.55,  1.3,  1.35 

[5] 

a,  ava2, 

and  a3 

Rate  of  moving  from  Hi,  H11,  Hie  and  Hir 
classes  to  H2,  H21 ,  H2C  and  H2R  classes 

0.0303,0.04,  0.16, 
0.05,  0.06 

[5] 

The  total  population  N  (t)  can  be  obtained  from 

m  =  S  (0  +  V  (0  +  Ht  (?)  +  H2(t)  +  Hu  (?)  +  H2I  (?)  +  Hlc  (?)  +  H2C  (?)  +  HlR  (?)  +  H2R  (?) 
+  IB(t)  +  CB(t)  +  RB(t) 


(14) 


n 


Here,  it  is  important  to  note  that  in  the  absence  of  the  disease  N(t)  — >  —  Moreover,  under  the 

Mo 

dynamics 

described  by  the  above  systems  of  equations,  the  region 

:  =  (S(t),V(t),Hl(t),H2(t).H„(t).H2,(‘).H,c(t).H,c(t),Ha(t),H2K(t),Il,(t),C,(t),N)eil'’\ 


Q  = 


S>  0,V>0,  Hx  >0,H2>  0,  Hu  >  0,H2I  >  0,  Hlc  >  0,  H2C>0,HlR>0,  H2R>  0,  IB>0,CB>0,N< 


jU  +  7T 

Mo 


(15) 

Is  positively  invariant.  Hence  the  system  is  both  mathematically  and  epidemiologically  well- 
posed.  Therefore,  for  initial  starting  point  x  e  id13 ,  the  trajectory  lies  in  Q.  Thus  we  restrict  our 
analysis  to  the  region  Q..  (Where  the  models  make  biological  sense) 
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Lemma  1  All  the  solution  of  the  Equations  2.1  -2.13  are  positive  for  all  time  t  >  0  provided  the 
initial  condition  are  positive. 

Proof:  Let 


(S(0)>V(0),H1(0),H2(0),/f„(0),H2,(0),/f1(;(0),ff!c(0),H1,(0),ffM(0),/1,(0),C1,(0),  AT)e< 


—  =  x  +  Ha>(l-vCt-vlHtc-v1Hx)+TV-(l-S)(l„+Xt)S-(v„+S)S 

>-(\-e)(x„+AB)s-(f,a+s)s 

This  implies, 

s'(,)>-(i-s)(ah+zb)s-(h,+s)s 


integrating  we  have 
similarly  it  can  be  shown  that 

S  >  0,  V  >  0,  Hl  >  0,  H2  >  0,  Hu  >  0,  H2I  >  0,  HIC  >  0,  H1C  >  0,  HlR  >  0,  H2R  >  0 ,IB>  0,  CB  >  0, 


(16) 

(17) 

(18) 


for  all  time  t  >  0.  Hence  all  solutions  of  the  HBV  model  remain  positive  for  all  non-negative  initial 
conditions. 


3.  Optimal  Control  Modeling 

ini 

—  =  7t  +  ju(d(1-v(1-u0)Cb-v1(1-u1)Hic-v2(1-u2)H2C)  +  tV-(1-S)(Ah+Ab)S-(ju0  +  S)S(19) 


^2 —  —  //(l -co)  +  SS  +  cjIb  —{Ah  +7tbAb)V-(ju0  +  r)y 

HTJ 

dH , 


dt 

dH^ 
dt 
dH \ 


dt 
dH , 


ccH i  (p2XBH2  (//0  T d H 2  , 

-  —  (l  —  o:)(l  —  w3 ) d-g/Zj  4- 1// 1  (l  —  Yy ) d-H/g  +  (f)\  (4  —  ccl ))HU  , 

-  —  (p2ABH2  +  cclHl  j  —(/j0  +  dlH  +  (j)2 )  H2J  , 


1  ,c 


dt 

dH2C 

dt 

dHlR 


dt 


dH 


2  R 


dt 

dlB 

dt 

dCR 

dt 


iU(0\'l  (1  )  Hl  c  +  (4  ctl )  Hl j  +  y/2p2AHCB  +  ^4// +  6^2  6r2 )) di\c  > 

=  //<»v2  (l  -  m2  )  H2C  +  q2(j)2H2I  +  a2Hlc  -  (//0  +  <Z2W  +  <92 )  H2C  , 

—  (l  _  ^1  )  (l  _  ^1  )  ^4./  “I"  ^1  (l  —  6C2  )  ^l,c  —  ( Ao  )  H \K  ■ 


^2  (l  ^2  )  ^2/  ^sHir  ( //Q  +  d3H  )  Z/2g  > 


=  (^ +(i-^)^)4-^(1-ri);iH4-(Ao+^+n)4 

■  —  pcov (4 —  u0 ) CB  +  q/il B  —ys 2  (l —  7*2 ) Cb  —  ( Ao  dB  +  ) CB 


;20) 

2i) 

;22) 

23) 

24) 

25) 

26) 

27) 

28) 

29) 

30) 
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dRB 

dt 


(l  q}yxIB  +  YtPb  Mo^b 


(31) 


Here,  the  control  parameters  employed  in  the  above  model;  Uo,  Ui,  U2  and  U3  are  the  effective  rate 
of  professional  effort  that  the  medical  expatriate  put  in  place  in  prevention  or  eradicate  of  mother 
to  child  transmission  in  the  HBV  carrier  class,  dually  infected  (HBY-HIY)  asymptomatic  carrier 
class,  dually  infected  (HBV-HIV)  symptomatic  carrier  class  respectively  while  U3  is  the  treatment 
rate  of  Hi,  1  dually  infected  (HBV/HIV)  asymptomatic  acute  class. 

Where 

Z0  —  Mo  +  S,  Zi  —  M 0  +T,  z2  —  Mo  +  CC,  Z2  —  Mo  dH ?  £4  —  Mo  (l  —  ^1 ) »  Z5  —  Mo  ^\h  +  » 

—  ^1^1  (l  —  ^1  )  >  Zj  ^8  —  Mo  +  ^4 II  +  ^1  (l  _  ^2  )  ’  ^8  —  Mo  ^2H  @ 2  ’ 

Zio  =^l(1-«l)(l-^l),  Zu  =6,l(1-«2)’  Zl2  =  M>+«3>  Zl3  =^(1-^).  Zl4  =A)+6?3H’ 

Zl5  =^(1-/l)’  *16  =  M)+0’  +  rp  *17  =^2(1-72)’  *18  =Mo+dB+r2- 

The  objective  function; 

Ch  B  B 

j (w0,w19w2,w3  ) — J  a^ic  (0 c  (0 (0 A^i/  (0 ~^~uo  (0 (0 

*^0  Z  Z 


+  ~^ul  (t)  +  ~^u2  (0  dt 

Using  similar  control  variables  and  objective  function  (as  in  above)  we  find  that  the  Hamiltonian 
takes  the  form; 


BBS 

H  =  AxHlc  (f)  +  AJi2C  (/)  + (f )  H-  (t^  +—u,Q  (t)  +—ux  (t)  +  —u2 
H=AlHlc(t)  +  A,Hlc(t)  +  A,C,(l)  +  A,Hll(,)  +  ^ul(t)  +  ^-ui(t)  +  ^-ul 


T>  10 

{t)  +  ^rul{t)  +  'YJJiifl 

^  i= 1 

(')+y«;(') 


(32) 


(33) 


+A4^ocHl  <p2ABH2  (//0  +  dH^H2~^+  X5  (l  ^)(l  «,)V^i+^i(l  Y\)^h^b  {Mo  +  <4  (l  ^i))^i/J 

MO)vi  (l —  Mj )  //j  c  +  q^Pi  (l  —  cq )  //[  7  +  y/2p2XHCB 

-  (mo  +  ^4//  +  a2  +  ^1  (1  "  ai  ))  Hlc 

\nBv+{\-5)s)xB-¥{\-rl)KiB 

-{Mo+^  +  h)h 

+Ao  |^/^®'/'(l  —  ^0  )  Cg  +  QYJb  ~V^2  (l  ~  Yl)  Cfl  —  (Mo  d  b  +  ^2  )  Cg  J 
The  adjoint  equations; 

A  =  +[(A  +^A)+Zi  ]4- Ai 


+  /lg|^/t/<uv,2  ( 1 —  u2  )  H2C  +  q2(p2H2I  +  a2Hlc  —  (  Mo  ^2 h  +  ^2 ) 


+  Z6[^2/LB//2  +0ri^l,/  “(Ml  +<^m  +^2)^2,/]+Z7 
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A/  =[^1(l-a)(l-M3)(A3  +Z2A3  -<xA4  -^(l-o^A^Aj 

A  =[%4+z3R-MA 

A5  =  -A4  +  Z4A5  -  ttjAg  +  ZgAy  +  ZjqA, 

A  =z^K~(iAK 

A/  =  -A,  +  //(wvj  (1-Mj)(Aj  -A7)  +  Z8A7  -a2Ag  -ZjjAg 

Ag  —  —  Aj  +  fj,(ov 2  (l  —  w2  )(Aj  —  Ag )  +  Z9Ag  —  ^2Aj0 

A,'  =  -crA, -^,(l-r,)AwA5+  (Z15  +ZI6)A1I  -37, Aj2 

A[0  =  — ^3  +  //®v(l  —  m0)(Aj  —  Aj2)  +  Z7AhA7  +Z17AhAj2  —  Z18Aj2 


With  transversality  conditions  A(  )  =  0  (i  =  1, 10,  tf  is  the  end  time).  Using  the 
Hamiltonian,  we  obtain  the  optimality  conditions; 

_  ('^l  _  ^12 

A 

(A1-A7)^v1i/1, 


M, 


U, 


11C 


W: 


_  (4  -^s)  7^®^ 2^ 2C 

=  ^ 

_  (A3  -^5)^1  (l-or)ABHi 


(37) 

(38) 

(39) 

(40) 

(41) 

(42) 

(43) 

(44) 


(45) 

(46) 

(47) 

(48) 


4.  Results  and  Discussion 

To  obtain  further  insight  into  the  influence  of  the  control  parameters  on  the  behavior  of  the 
systems  under  consideration,  thus,  we  consider  the  following  Figures  plotted  with  the  aid  of  Maple 
software  [7]. 


- Asymptomatic  HIY.  Hl(t)  (without  control) 

- Symptomatic  HIV  H2(t)  (without  control) 


- .Acute  HBV  (without  control) 

Carrier  HBV (without  control) 


Figure  2 


Figure  3 
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- I_B(t) - C_B(t)  H_l(t) - H_2  (t)  | 

Figure  4 


- H_lC(t)  showing  (u_3=0,  0.3,  0.6,  0.9) 

- H_1 1(t)  showing  (u_3=0,  0.3,  0.6,  0.9) 

H_2C(t)  showing  (u_3=0,  0.3,  0.6,  0.9) 
H_2I(t)  showing  (u_3=0,  0.3,  0.6,  0.9) 


Figure  5 


Figure  6  Figure  7 


- R_B(t)  show  ing  (ul,  u2,  n3=0,  0.3,  0,6,  0.9) 

- HlR(t)  showing  (ul,  u2,  u3=0,  0.3,  0,6,  0.9 

H2R(t)  showing  (ul,  u2,  u3=0,  0.3,  0,6,  0.9) 


Figure  8 
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In  the  co-infection  model  above,  the  control  parameters  employed  are;  Uo,  Ui,  U2  and  U3  are  the 
effective  rate  of  professional  effort  that  the  medical  expatriate  put  in  place  in  prevention  or 
eradicate  of  mother  to  child  transmission  in  the  HBV  carrier  class,  dually  infected  (HBV-HIV) 
asymptomatic  carrier  class,  dually  infected  (HBY-HIV)  symptomatic  carrier  class  respectively. 
While  U3  is  the  treatment  rate  of  HIV  asymptomatic  acute  class  (Hi)  without  considering  any 
optimal  control  measure  except  the  therapeutic  vaccination.  We  considered  the  state  of 
asymptomatic  HIV  (Hi)  and  symptomatic  HIV  (H2)  (without  control)  in  Figure  2.  Similarly,  in 
Figure  3,  relationship  between  acute  infection,  HBV  (Ib)  and  Carrier  HBV  (Cb)  were  considered. 
Figures  4,  harbored  Figures  2  and  3.  Increasing  the  controls  (Figures  5)  from  0,  30,  60  and  90 
percents,  reduces  the  dual  infection  in  Hie,  H11  and  H2C  respectively. 

Thus,  Hie  show  major  significant  difference  in  response  to  the  control  adopted,  while  H21  made  no 
difference  and  die-out.  Figures  6  shows  the  response  (decrease)  of  Hie  and  H2C  with  the  increase 
in  control  measure.  Also,  Hie  registered  a  major  significant  difference  when  compared  with  H2C. 
Meanwhile,  Figure  7  displays  the  relationship  between  the  Susceptible  and  Vaccinated  class, 
which  decreases  and  moves  at  a  constant  rate.  In  addition,  the  sensitivities  of  the  recovery  classes 
are  obtained  in  Figure  8,  where  Hir  recorded  a  more  pronounced  significant  difference,  also  in 
H2R,  however,  little  or  none  is  noticed  with  Rb.  A  similar  pattern  of  some  of  the  results  was 
obtained  in  [6],  however,  when  comparing  with  other  results  in  literature  we  are  able  to  record  a 
significant  changes  in  the  behavioral  dynamics  of  the  model  owing  to  the  sensitive  control 
measure  put  in  place. 

5.  Conclusion 

This  scheme  has  shown  that  the  disease  can  be  resisted  with  vaccination(s)  and  treatment.  As  the 
results  mentioned  earlier  reveal,  with  a  two-year  robust  vaccination  campaign,  the  epidemic  will 
be  avoided,  but  the  rate  drops  only  to  move  at  a  constant  rate.  However,  we  must  be  cognizant  of 
the  fact  that  a  two-year  campaign  does  not  only  eliminate  future  epidemics,  but  the  routine  must 
be  repeated  to  avoid  future  occurrence.  So  the  more  active  our  campaign  strategies  (vaccination(s) 
and  treatments)  are,  the  sooner  the  future  epidemic  will  be  arrested. 
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